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We study the two point correlation function of two nucleon currents 
sandwiched between the vacuum and the rho meson state. The 
light cone QCD sum rules are derived for the pNN vector and 
tensor couplings simultaneously. The contribution from the excited 
states and the continuum is subtracted cleanly through the double 
Borel transform with respect to the two external momenta, pf, p\ = 
(SI ; (p — q) 2 . Our results are: g p = 2.5 ± 0.2, n p = (8.0 ± 2.0), in good 

agreement with the values used in the nuclear forces. 
PACS Indices: 21.30. +y; 13.75. Cs; 12.38.Lg 
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1 Introduction 



Quantum Chromodynamics (QCD) is asymptotically free and its high energy behavior 
has been tested to two-loop accuracy. On the other hand, the low-energy behavior has 
become a very active research field in the past years. Various hadronic resonances act as 
suitable labs for exploring the nonperturbative QCD. Among which, the inner structure of 
[ nucleon and mesons and their interactions is of central importance in nuclear and particle 
c3 physics. 

Internationally there are a number of experimental collaborations, like T JNAL (former 
CEBAF), COSY, ELSA (Bonn), MAMI (Mainz) and Spring8 (Japan), focusing on the 
nonperturbative QCD dynamics. Especially the Mainz research project MAMI (Mainzer 
Mikrotron) with its planned extension from a 855 MeV to 1.5 GeV electron c.w. acceler- 
ator and the Japan Hadron Facility (JHF) at Spring8 will extensively study the photo- 
and electro-production of vector mesons off nucleons. 

Moreover, the strong pNN couplings are , like ttNN and ttNA couplings, the basic 
inputs for the description of nuclear forces in terms of meson exchange between nucleons. 
So far the linkage between the underlying theory QCD and the phenomenological pNN 
couplings has not been made. Especially the commonly adopted tensor-vector ratio k p = 
6.8 is much larger than the vector meson dominance model (VDM) result k v = 3.7. 
One wonders whether it is feasible to calculate the pNN couplings directly with the 
fundamental theory QCD? 
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Although it is widely accepted that QCD is the underlying theory of the strong in- 
teraction, the self-interaction of the gluons causes the infrared behavior and the vacuum 
of QCD highly nontrivial. In the typical hadronic scale QCD is nonperturbative which 
makes the first principle calculation of these couplings unrealistic except the lattice QCD 
approach, which is very computer time consuming. Therefore, a quantitative calculation 
of the pNN couplings with a tractable and reliable theoretical approach proves valuable. 

The method of QCD sum rules (QSR), as proposed originally by Shifman, Vainshtein, 
and Zakharov and adopted, or extended, by many others 0, ^ ffl, are very useful in 
extracting the low-lying hadron masses and couplings. In the QCD sum rule approach the 
nonperturbative QCD effects are taken into account through various condensates in the 
nontrivial QCD vacuum. A recent review of QSR is given by Shifman [||. In this work 
we shall use the light cone QCD sum rules (LCQSR) to calculate the pNN couplings. 

The LCQSR is quite different from the conventional QSR, which is based on the 
short-distance operator product expansion. The LCQSR is based on the OPE on the 
light cone, which is the expansion over the twists of the operators. The main contribution 
comes from the lowest twist operator. Matrix elements of nonlocal operators sandwiched 
between a hadronic state and the vacuum defines the hadron wave functions. When the 
LCQSR is used to calculate the coupling constant, the double Borel transformation is 
always invoked so that the excited states and the continuum contribution can be treated 
quite nicely. Moreover, the final sum rule depends only on the value of the hadron wave 
function at a specific point, which is much better known than the whole wave function H]. 
In the present case our sum rules involve with the rho wave function (RWF) (p p (u = |). 
Reviews of the method of LCQSR can be found in 0, §. 

The LCQSR has been widely used to treat the couplings of pions with hadrons. 
Recently the couplings of pions with heavy mesons in full QCD ||, in the limit of tuq — > oo 
P], 1/rriQ corrections and mixing effects |]10|| , the couplings of pions with heavy baryons 
jrj, the ttNN and 7riViV*(1535) couplings the p -> tttt and K* -> Kir decays 
and various semileptonic decays of heavy mesons are discussed. 

The QCD sum rules were used to analyze the exclusive radiative .B-decays with the 
help of the light-cone vector meson wave function in ||15|| . With the same formalism the 
off-shell Qb*b p and gn*D P couplings in fl6| and the p decay widths of excited heavy mesons 
| 17| were calculated. 

Our paper is organized as follows. Section [l] is an introduction. We introduce the 
two point function for the pNN vertex and saturate it with nucleon intermediate states 
in section ||. The definitions of the RWFs and the formalism of LCQSR are presented in 
section [3[ In the next section we present the LCQSR for the pNN coupling. In section 
|5] we present some discussions of these RWFs and their values at the point u = ~. We 
make the numerical analysis and a short discussion in section || 



2 Two Point Correlation Function for the pNN coupling 

Many authors have studied the strong pNN couplings. It was pointed out that the 
inclusion of an effective p-pole contribution leads to a large value for the tensor-vector 
coupling ratio k p = 6.6 ± 1 JTB], [T!|, in the dispersion-theoretical analysis of the nucleon 
electromagnetic form factors. The above value is consistent with other determination 



20| , plf . Brown and Machleidt have discussed the evidence for a strong pNN coupling 
from the measurement of e\ parameter in NN scattering Brown, Rho and Weise 

suggested that consistent with a quark core radius of 0.5fm for the nucleon 



2 



and an equal factorization of the baryon charge between the quark and meson cloud in the 
two-phase Skyrme model p3| . Recently Wen and Hwang used the external field method in 
QCD sum rules to study the pNN couplings |23|]. They obtained k p = 3.6, in agreement 
with VDM result k v = 3.7 and k s = —0.12. In their work the authors introduced the 
vector-like p-quark interaction Lagrangian by hand and treated the vector meson quark 
coupling as free parameter. In other words, the pNN vector and tensor couplings can not 
be determined simultaneously. 

We shall calculate the pNN vector and tensor couplings simultaneously using vector 
meson light cone wave functions up to twist four, which will result in a reliable extraction 
of K p . 

We start with the two point function 

n(pi,p 2 ,g) = / d*xe^(0\T Vp (x)ff p (0)\p O (q)) (1) 
with pi = p, p 2 = p — q and the Ioffe nucleon interpolating field || 

V P (x) = e abc [u a (x)C lp u b (x)] l5 rd c (x) , (2) 

fj P (y) = e abc [u\y) lu Cu aT (y)]d c (y)Yl* , (3) 

where a,b,c is the color indices and C = 27270 is the charge conjugation matrix. For the 
neutron interpolating field, u «-> d. 

The rho nucleon couplings are defined by the pNN interaction: 

£ p nn = g p( fN[i, + Kp^W ■ (4) 

p p is an isovetor in (^). g p is the rho- nucleon vector coupling constant and k p is the 
tensor-vector ratio. 

At the phenomenological level the eq.(|l|) can be expressed as: 

U(p uP2 , q) = i\ N g p e^(X) (pl-m%)te-m%) + ' ' ' (5) 

where e M (A) is the rho meson polarization vector. The ellipse denotes the continuum and 
the single pole excited states to nucleon transition contribution. is the overlapping 
amplitude of the interpolating current r] N (x) with the nucleon state 

(0\ V (0)\N(p)) = X N u N (p) (6) 

Expanding @ with the independent variables P = pi ^ P2 , q and decomposing it into 
the chiral odd and chiral even part, we arrive at: 

n = n D + n e , (7) 

where 

U ^ q) = (pl-^- m ^ • P)P + ^ 

_ l + 2^ (e _ ^ + . (l + e spagP.fy.y + . . . j (8) 
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and 



n«(P, P2 ,,) = (ri - m Srf-,»&) {(2ro " + ir/^ ■ p ) 

1 + ftp /AA „„ s K 



m N (eq - qe) - - (e • P)(qp - pq)} H , (9) 

with q 2 = m 2 . We have not kept the single pole terms in (H) and (§) since they are always 
eliminated after making double Borel transformation in deriving final LCQSRs. 

It was well known that the sum rules derived from the chiral odd tensor structure yield 
better results than those from the chiral even ones in the QSR analysis of the nucleon 
mass and magnetic moment [|3], We shall consider chiral odd tensor structures only 
below. 



3 The Formalism of LCQSR and Rho Wave Functions 

Neglecting the four particle component of the rho wave function, the expression for 
U(p1,P2, q 2 ) with the tensor structure at the quark level reads, 

J e^dx(0\T Vp (x)f} p (0)\p°(q)) = -4e abc e a ' b ' c ' W S a /(x) 

7 ^(0| C / c (x)« c '(0)|p + (g)) T C 7M ^r'(x)7,75 , (10) 

where iS(x) is the full light quark propagator with both perturbative term and contribu- 
tion from vacuum fields. 

iS(x) = (0\T[q(x), q(0)]\0) = - M _ ^-(qg s a • Gq) 

+ ^Sh t du{2(l-2u)x^ u + ie IMUpal5 Yx a }G> lu + ■■■ . (11) 

lD7T z JO 

where we have introduced x = £p7 M . In our calculation we take the tiny current quark 
mass to be zero. 

The relevant feynman diagrams are presented in FIG 1. The squares denote the rho 
wave function (RWF). The broken solid line, broken curly line and a broken solid line with 
a curly line attached in the middle stands for the quark condensate, gluon condensate and 
quark gluon mixed condensate respectively. 

By the operator expansion on the light-cone the matrix element of the nonlocal op- 
erators between the vacuum and rho meson defines the two particle rho wave function. 
Up to twist four the Dirac components of this wave function can be written as follows 
. For the longitudinally polarized rho mesons, 

p(A) . T »i m 2 x 2 
<0|fi(0)7^(x)|p-(ff,A)> = / p m / ,{( 6i — - / due- lu ^(u^ 2 ) + ^—A(u^)} 

C£ X J 

+ (4 A) - q, ) / due-^ x g\u^ 2 ) - -x,- - 2 m 2 p / du e~ m ^C{u, /i 2 )} , (12) 

p qx Jo 2 [q ■ x) A H Jo 

and 

(0\u(0) lfll5 d(x)\p-(q,X)) = -lf p m p e^ef v q a x p £ due^g^u, // 2 ) (13) 

where 

C{u)=g,{u) + ^{u)-2g { l\u). (14) 
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The link operators Pexp ig J dax^A^ax) are understood in between the quark 
fields. The distribution amplitudes describe the probability amplitudes to find the p in 
a state with quark and antiquark carrying momentum fractions u (quark) and 1 — u 
(antiquark), respectively, and have a small transverse separation of order - pTfl . 

The vector and tensor decay constants f p and fj are defined as 

(0|u(0) 7 ^(0)|p-(g,A)) = f p m p ef . (15) 
All distributions <fi = {<p\\,g v ,g a ,A,C} are normalized as 

du(f)(u) = l, (16) 



The twist-three three-particle quark-antiquark-gluon distributions are | 27| : 

(0\u{0)-f a g s G pu {ux)d(x)\p-(q,\)) = iq a [q p e^ } u - q v e%]f^ p V{u,qz) + ■■■ , (17) 

(0\u(0)^ 5 g s G l , u (ux)d(x)\p~(q, A)) = q a [q v eY p - q^l\f 3p A{u, qx) + ■ ■ ■ , (18) 

where the operator G a/ 3 is the dual of G a p\ G a/ 3 = \e af 3s p G 5p and the ellipses denote 
higher twists contribution. 

The following shorthand notation for the integrals defining three-particle distribution 
amplitudes is used: 

T(u,qx)= fvae- iq < a3+ua ^F(ai,a g ,a 3 ). (19) 

Here T = {V, ^4} refers to the vector and axial- vector distributions, a is the set of three 
momentum fractions: «3 (d quark), ot\ (u quark) and a g (gluon), and the integration 
measure is defined as 

J Va = J dai J da 3 J da q 5(1 — ^2 a i)- (20) 

The normalization constants f^ p , f^ p , are defined in such a way that 

fVa(a 3 -a 1 )V(a 1 ,a g ,a 3 ) = 1, (21) 

JvaA(a u ag,a 3 ) = 1. (22) 

The function A is symmetric and the functions V is antisymmetric under the interchange 
a\ <-> «3 in the SU(3) limit (|28|, |2~7fl), which follows from the G-parity transformation 
property of the corresponding matrix elements. 

In the infinite momentum frame the RWFs are associated with the leading twist 
two operator, A(u) correspond to twist four operators, and g a (u), g v {u) to twist three 
ones. The three particle RWFs V, A are of twist three. Details can be found in [p7| . 



4 The LCQSR for the pNN coupling 

Expressing ( |10|) with the longitudinal rho wave functions (LRWFs), we can obtain the 
expressions for the correlator in the coordinate space. Up to dimension six the gluon 
condensate is the only relevant condensate contributing to the chiral odd tensor structure. 
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Yet the gluon condensate always appears with a large suppression factor, which arises from 
the two- loop internal momentum integration in the diagram (d) in FIG 1. Its contribution 
is quite small, which is confirmed by our detailed calculation. For example, after double 
Borel transformation, diagram (d) is suppress by a factor where (g 2 s G 2 ) = 0.48GeV 4 

and Ml ~ 1.5GeV 2 . 

Diagram (a) involves with two-particle LRWFs. After tedious but straightforward 
calculation, we get: 

Il2(pi,P2,g) = / d A x £ due^- u ^f p m p {-£^e aPlMT e^q^j a 75 



tv 4 x 6 (q ■ x) 7r 4 x 6 p (q ■ x) 



-Al'»;*»«R - ^MG')^"t^ + ■■•>. (23) 

The LRWFs can be found in the previous section. Diagram (b) is associated with vac- 
uum gluon fields. But its contribution vanishes due to isospin symmetry. Our explicit 
calculation confirms it. 

We frequently use integration by parts to absorb the factors l/(q ■ x) and l/(q ■ x) 2 , 
which leads to the integration of RWFs. For example, 

/ Tp(u)du = i / e- luq - x ^(u)du + ^(u)e- luqx \ l Q , (24) 

Jo q ■ x J 

where the functions ty(u) is defined as: 

\&(u) = + f U ip(u)du . (25) 
Jo 

Note the second term in (H) vanishes after double Borel transformation or due to 4> p {uq) = 
^(uo) = at end points uq = 0,1. 

We first finish Fourier transformation. The formulas are: 

-d D x -. 2D ~2m!L F *' P / 2 7 — ■ (26) 



{x 2 ) n V ' (-pl)D/2-n p( n ) 

xe ipx D , n+1 2 D - 2n+1 p°/ 2 T(D/2 + 1 
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X 



2\n 



ft-(-'r (Y) „ ; ;,., ' rw ' *< (27 » 



x p x u e ipx D 1 \n p-2n+i p/2r 9fw T(D/2 + l n) 

{X 2 Y l(_p2)U/2+l-n r(n ) 

2p pVv r(D/2 + 2-n) 

+ (_ p 2)D/2+2-n Y {n) 1 ' 1 ' 

The next step is to make double Borel transformation with the variables p\ and p 2 to 
(|J) and (]30|)-(|3T|). The single-pole terms in (|5]) are eliminated. The formula reads: 

B x M h2 M K , 1 sr = (M 2 ) 2 -e-^5(n - u ) , (29) 

hi 2 v\ [ m 2 _ (1 _ u )p\ - upl] n y ' K ' ' V ; 

where w = M i^ M i , ^f 2 = yfc+M* ' M% are the Borel parameters. 

Finally we identify the same tensor structures both at the hadronic level and the quark 
gluon level. Subtracting the continuum contribution which is modeled by the dispersion 
integral in the region s\, S2 > Sq, we arrive at: 







ie 

„2 



W2c/ P — ^m p e ^ = -|-^e m 2 {g {u )M M^. 



-mjG 3 MMVi(^) + ^(^G 2 )M 2 / o(^)} , (30) 



i(e-P)P 



where g 2 = m 2 , = 1 — e _a: |j is the factor used to subtract the continuum, s 

k=0 

is the continuum threshold. The sum rules are symmetric with the Borel parameters M 2 
and Mf. It's natural to adopt Mf = Aff = 2M 2 , m = §. The functions Gi(u),i = 
0, 1, 2, 3, are defined as: 



ru 

G {u) = J dt(f>\\{t) , 


(31) 


G 1 {u) = f U dtg v (t), 
Jo 


(32) 


ru 

G 2 (u) = / dtA(t) , 
Jo 


(33) 


ru rt 

(it) = dt dsC(s) . 


(34) 


Jo Jo 


and parameters 





The resulting sum rules depend on the RWFs, the integrals and derivatives of them at 
the point Uq = |. The distribution amplitudes of vector mesons have been studied in 
33, Ea ETfl using QCD sum rules. We adopt the model RWFs in Refs.pl. 



The values of the two-particle RWFs, their derivatives and integrals at the point 
M = \ using the form in |9[ are: g a = 1.15 ±0.23, g v = 0.64, (f>\\(u ) = 1.1, G (u Q ) = 0.5, 
Gi(u ) = 0.5, G 2 (u ) = 0.58, G 3 {u ) = -0.13, A(u ) = 2.18. 

The experimental value for the rho meson mass m p and the decay constant f p is 
f p = 198 ± 7MeV and m p = 770MeV || . 

The various parameters which we adopt are so = 2.25GeV 2 , rriN = 0.938GeV, Xn = 
0.026GeV 3 at the scale \x = lGeV. The working interval for analyzing the QCD sum 
rules for nucleons is 0.9GeV 2 < M§ < 1.8GeV 2 , a standard choice for analyzing the 
various QCD sum rules associated with the nucleon. 



6 Numerical analysis and results 

In order to diminish the uncertainty due to A^v, we shall divide our sum rules by the 
famous Ioffe's mass sum rule for the nucleon: 

32n*\%e~& = M 6 f 2 A + -M 2 / ( *L) + t a * _ . (35 ) 

n J y M 2 ' A y M 2 ' 3 3M 2 v ; 

Dividing (]3Qf)-d3S|) by (|3~5"1), we have two new sum rules for g p (l + K p ),g p . The de- 
pendence on the continuum threshold Sq and Borel parameter M 2 of these sum rules are 
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presented in FIG 2-3. From top to bottom the curves correspond to Sq = 2.35,2.25,2.15 
respectively. These sum rules are stable with reasonable variations of So and M 2 as can 
be seen in FIG 2-3. Numerically we have: 



for 



for || 



and 



g p {l + n p ) = (22 ± 3) , (36) 
g p = (2.5 ± 0.2) , (37) 



We can also divide ( j30|) by ([H]) to get a new sum rule for 1 + n p . 

P ~ mJGs(«o)/o(^r) ' [ } 

The result is presented in FIG 4. Numerically, 

1 + k p = (9.0 ± 2.0) , (39) 

which corresponds to 

k p = (8.0 ± 2.0) . (40) 
Brown and Machleidt emphasized that the strong pNN coupling 

g£l±^! = (37 ±13) (41) 

47T 

k p = (6.6±0.1) (42) 



should be adopted in order to reproduce experimental data [p^] . The vector meson dom- 
inance (VMD) model yields 



Air 

Our result is 



13.25 . (43) 



,_ =(39 ±10), (44) 

which agrees very well with (|41|) and deviates strongly from VMD prediction (f43|). 

We have included the uncertainty due to the variation of the continuum threshold 
and the Borel parameter M 2 in our analysis. Other sources of uncertainty include: (1) 
the truncation of OPE on the light cone and keeping only the few lowest twist operators; 
(2) the inherent uncertainty due to the model RWFs etc. In the present case the major 
uncertainty comes from the RWFs since our final sum rules depends both on the value of 
RWFs and their integrals at u . 



Our result k p = (8.0±2.0) is about two times larger than that derived in pi| , k p = 3.6 
treating the rho meson as the external field. The reason is two-old. Firstly, the vector- 



like pqq interaction is assumed in [23|. Moreover, for the antisymmetric sum rules the 
susceptibilities Xi K an d £, take the same values as in the nucleon magnetic sum rules 
where the electromagnetic field is treated as the external field. In our opinion, such a 
treatment employs the vector meson dominance assumption inexplicitly, which may be 

the underlying reason for Kip ~ Kjy 3.7. Secondly, the vector meson mass corrections 

turn out to be large as pointed out in |p9[ . This effect is explicitly taken into account in 
our calculation. 
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In summary we have calculated the pNN couplings starting from QCD. The contin- 
uum and the excited states contribution is subtracted rather cleanly through the double 
Borel transformation in both cases. Our result strongly supports large value for the 
tensor-vector ratio k p in the nuclear force. 

Acknowledgments: This work was supported by the Postdoctoral Science Foundation 
of China and Natural Science Foundation of China. 
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Figure Captions 



FIG 1. The relevant feynman diagrams for the derivation of the LCQSR for pNN 
coupling. The squares denote the rho wave function (RWF). The broken solid line, broken 
curly line and a broken solid line with a curly line attached in the middle stands for the 
quark condensate, gluon condensate and quark gluon mixed condensate respectively. 

FIG 2. The sum rule for g p (l + a function of the Borel parameter M 2 for ** 

with the model RWFs in pg[ |. From bottom to top the curves correspond to the continuum 
threshold s = 2.35, 2.25, 2.15GeV 2 . 

FIG 3. The sum rule for g p as a function of M 2 and so for 
FIG 4. The sum rule for 1 + k p as a function of M 2 and so from 
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